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Abstract. We have embarked upon a project to model the UV spectra of BALQSOs 
using a Monte Carlo radiative transfer code p reviously validated th rough modelling of 
the winds of cataclysmic variable stars (e.g. iNoebauer et aDl2010h . We intend to use 
the simulations to investigate the plausibility of geometric unification (e.g. lElvisll2000l) 
of the different classes of QSO. Here we introduce the code, and present some initial 
results. These demonstrate that for reasonable geometries and mass loss rates we are 
able to produce synthetic spectra which reproduce the important features of observed 
BALQSO spectra. 



1. Motivation 

Broad absorption line quasars (BALQSOs) are a class of QSO exhibiting broad, blue 
shifted absorption lines of highly ionized species such as N v, C iv, O vi and Si iv . 
These objects make up perhaps 17% of the population of QSOs (e.g. lKnigge et al.ll2008l : 



iHewett & Foltdl2003h . Their unusual UV spectra are be lieved to be the res ult of fast 
(0.01-0.2c) winds originating from the accretion disk (e.g. lKorista et al.|[T9 92). although 



the geometry and kinematics of these winds are not well understood, and even the 
driving mechanism is unclear. It is a source of discussion as to whether this observation 
demonstrates that all QSOs spend 17 % of their lifetime as a BALQSO ('evolutionary 



unification'; e.g. [Becker etaDl2000h . or whether all QSOs have these observational 



sign atures when observed from 17% of possible directions ('geometrical unification'; 
e.g. lElvisll2000b . Our aim in developing a code to simulate UV spectra for quasar disk 



winds is to permit quantitative evaluation of prospective geometries for geometrical 
unification. 



2. Python - A Radiative Transfer code 



Python is a Monte Carlo type radiative transfer code, first described in lLong & Knigge 
(2002 ), which has been use d previously to simulate the spectra of cataclysmic v ariables 
(e.g. INoebauer et ai1l2010h and massive young stellar objects dSim et al.ll2005h . Cata- 
clysmic variables are thought to drive winds with similar geometries, kinematics and 



1 



2 



Nick Higginbottom, Knox S. Long, Christian Knigge and Stuart A. Sim 



Initialization 



Get user inputs 



Generate initial 
temperature, density 
and velocity structure 
of the wind 



Ionization calculations 



3: 



Generate photon 
distributions 



Transport photons 
through wind 



T 



Compute temperature 
structure to balance 
I heating and cooling , 



Compute a new 
ionization state 



Spectrum simulation 

~ 1 



Generate photon 
distribution 



Transport photons 
through wind 



Calculate photon 
weight escaping in 
required observer 
directions 



Compute observed 
spectra 



Figure 1 . Schematic of Python execution cycle 



ionization conditions to those in AGN and are observed to produce similar wind-formed 
spectral lines and line profile shapes. It is therefore natural to attempt to apply Python 
to the modelling of BALQSO spectra and hence, if geometrical unification is possible, 
QSO spectra in general. 

Python is designed to be easy to adapt to different kinematic descriptions of winds. 
A number of azimuthally symmetric wind models are currentl y impl e mented in the 
code, including the wind prescriptions of IShlosman & Vitellol (Il993h . iKnigge et al.1 
(Il995h . and one simulating the wind envisioned by lElvisI (l2000h . 

Figure Q] shows schematically the execution stages of the code. Once the wind 
model is defined, and the computation grid populated by wind cells of given density 
and velocity, the code proceeds by generating a population of photons arising from all 
radiating parts of the simulated system. These can include: 

• the disk - a thin disk approximation is used here, implemented as a series of 
annuli radiating as black bodies (in principle any spectral energy distribution can 
be used); 

• the wind - radiation is produced via line and continuum processes at the current 
electron temperature; 

• a spherical 'corona' at the centre of the geometry - producing photons defined by 
a power law extending into the x-ray regime. In the current implementation this 
region is contained within the innermost extent of the disk. 

The initial ionization state of the wind is calculated using the Saha equation and an 
initial estimate of the electron temperature. Photons are then tracked through the wind. 
These photons are used to calculate the local radiation field for each cell, and this is then 
used to compute new ionization state. This process is repeated until the temperature and 
ionization state have converged throughout the entire outflow. 

The ionization calculations are optimised for speed to allow even complex wind 
models and radiation sources to be considered. Key to this is a parametrisation of 
the photon field in any given cell. This parametrisation is then used to solve for the 
ionization state via a modified version of the Saha equation. All of the previous work 
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using Python has been based on a dilute blackbody approximation (IMazzali & Lucy 
1993). However, it is unclear whether this method is appropriate for an AGN exhibiting 



a strong power-law component. Therefore we hav e recently incor porated a broken 
power-law approximation to the local radiation field (ISim et al.ll2008h into Python. We 
are currently experimenting with this to determine whether it will improve the estimate 
of the ionization state in AGN winds. 



3. Wind Geometry 

Here we present simulated AGN spectra based on the geometry suggested by lElvisI 
(l2000h . A schematic of the wind parameters is given in Figure[2|together with the values 
of the most important model parameters. Parameters controlling the radial dependance 
of the mass loading of the wind and the wind acceleration are omitted for brevity. 



Disk Parameters 


Central object parameters 

mass - 1O 8 M 
radius - 8.8 x 10 13 cm(6R G ) 
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Figure 2. Cartoon of the wind geometry for one possible Elvis-style wind 



4. Example results 

In Figure [3] we present some results from the geometry described above, illuminated 
with photons generated from the disk only. The top left plot shows the assigned wind 
velocity, and the the electron temperature calculated for the wind in the top right plot. 

The bottom left plot shows the calculated fractional abundance of C iv 1550A (the 
most ubiquitous and often strongest BAL feature) and the fracional abundance of O vi 
1240 A is given in the bottom right plot. These plots show that, in this model, O vi is 
found in the lower part of the wind, where the temperature is high, whilst we find C iv 
in the upper, cooler part of the wind. Also note that the velocity of the wind is higher 
in the upper part of the wind where the C iv is found. 

In Figure |4] a we present simulated spectra for observers with inclination angles 
i=0°, 63°, 69° and 85° for this wind model. 

The i=0° plot shows the calculated spectrum looking down on the disk/wind sys- 
tem, i.e. with the disk seen face on. The wind produces several emission lines, but 
the bulk of the emission is coming directly from the disk, the central parts of which 
are unobscured from this direction. For i=63°, the observers line of sight towards the 
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Figure 3. Wind velocity, calculated electron temperatures, C iv and O vi abun- 
dances for the model shown in Figure [2 

centre of the accretion disk lies entirely within the outflow. The wind accelerates along 
streamlines within this wind cone, and so material at a range of velocities is in the line 
of sight. C iv is found along much of the top edge of the wind, and so the 1550A line is 
visible at a range of velocities, giving rise to the characteristic deep, broad, absorption 
feature extending blue-ward of the C iv rest wavelength. As shown in Figured O vi is 
found in regions with lower velocity. The shape of the O vi 1240 A line reflects this: 
it is much narrower than the carbon line. Features due to other ions commonly used 
as BALQSO tracers are also seen with similar profiles. The dramatic effect on the line 
profiles of small changes in viewing angle is evident by comparing the i=63° spectrum 
with that computed for i=69°. For i=85° the observer is looking under the wind cone, 
and the C iv line, amongst others, is seen in emission. The line is broadened due to the 
Keplarian rotation of the wind close to the disk. 



5. Summary and future work 

We have begun a project to use Python, a Monte Carlo radiative transfer code, to simu- 
late theUVsrjectra of QSO winds. We present initial results for a geometry suggested 
by lElvisI (l2000h . illuminated by a disk only, and observe that the simulated spectra 
exhibit features reminiscent of those seen in BALQSOs. 

At present, for the model parameters we have tried, the wind becomes over-ionized 
whenever a strong central X-ray source is included, and none of the expected UV lines 
are observed. We are incorporating a dense shielding region into the geometry to see if 
this produces BAL-like spectra even with a strong central X-ray source. 
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Figure 4. Spectra for i=0°, 63°, 69° and 85° calculated for the model shown in 
figure [2] 



We will produce synthetic spectra for a range of geometries, viewed from different 
sight-lines to and so inform the discussion of the evolutionary vs orientation interpreta- 
tion of the BALQSO phenomenon. 
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